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Abstract 


Improved knowledge of the earth's gravity field can be obtained from new 
and improved satellite measurements such as satellite to satellite tracking 
and gradiometry. This improvement has been examined by estimating the 
accuracy of the determination of mean anomalies and mean undulations in 
various size blocks based on an assumed mission. In this report the accuracy 
is considered through a commission error due to measurement noise propa- 
gation and a truncation error due to unobservable higher degree terms in 
the geopotential. To do this the spectrum of the measurement has been related 
to the spectrum of the disturbing potential of the earth's gravity field. Equations 
were derived for a low-low (radial or horizontal seperation) mission and a grad- 
iomete r miss ion. 

For a low -low mission of six month's duration, at an altitude of 160 km, 
with a data noise of - 1 /im/sec for a four second integration time, we would 
expect to determine 1° x 1° mean anomalies to an accuracy of ±2.3 mgals and 
1° x 1° mean geoid undulations to ±4. 3 cm. A very fast Fortran program is 
available to study various mission configurations and block sizes. 
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Introduction 


A number of studies have recently been carried out to investigate the 
accuracy of the determination of the earth's gravity field from a gravity mapping 
satellite mission (or GRAVSAT). Such a mission has been discussed extensively 
for the past several years with a summary of applications and questions being 
found in the National Academy of Science report: Applications of a Dedicated 
Gravitational Satellite Mission. 

The studies performed to date include the following: Douglas et al. (1980) 
who describe an error analysis where mean gravity anomalies are unknown to 
be solved for in a least squares adjustment procedure; Rapp and Hajela (1979) 
who describe the use of least square collocation methods to a high-low satellite- 
to-satellite tracking mission; Breakwell (1980a, b) who describes an optimum 
filter approach in the frequency domain making a reasonable flat earth approxi- 
mation; Pisacane and Yionoulis (1980) who describe gravity improvement tests 
made with ORAN (Orbit Analysis Program) program; Lancaster et al. (1980) 
who describe error analysis studies for one degree mean anomalies using an 
adjustment procedure; and Rummel (1980) who describes the accuracies achievable 
from a low-low SST mission using least squares collocation methods. 

The above, and other, studies have given much additional insight into what 
might be obtained from a gravity field mapping mission. However, the studies 
still have not given a wholly satisfactory picture. For example, some studies 
only look at SST data; some solve for anomaly accuracies only; some use 
limited data sets; some use flat earth approximations. We wished to develop a 
rapid error analysis procedure that could be used for SST missions and gradi- 
o me ter missions; that could be used for mean anomaly and mean undulation 
accuracy estimation; and that could be used with almost continuous data. This 
led to the analysis and results described in this report. We do not believe this 
work cannot be improved upon; there will be areas where assumptions are made 
and missions not considered that could be (perhaps) worked with at a later time. 

We do feel the method presented here is reasonable and easy to use on computer 
implementation. 


The Method 


The basic method used here is the analysis of the gravity field mapping in 
the frequency domain with final results discussed in terms of mean gravity anomalies 
and mean geoid undulations on the surface of the earth. The types of measurements 
to be considered are velocity (or velocity differences) and vertical component 
gradiometer measurements. 
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The errors in the measurements are assumed to be uncorrelated. This 
type of random process can be approximated as white noise, although this is 
somewhat unrealistic since pure white noise is associated with infinite variance 
as the frequencies are affected equally throughout their infinite range. However, 
for the present purposes, it will be postulated that the error spectrum is flat at 
least over the frequencies of the measured data. 

To derive the error spectrum, we suppose that the physical quantity to be 
measured, say F, can be expressed at the satellite altitude as a series of 
spherical harmonic functions. Therefore, let 

00 

F( 0 , A) = E E fn B Y n8 (0,X) (1) 

n=0 n 

where (0,A) are spherical coordinates and the Y nm are fully normalized 
harmonic functions satisfying the orthogonality relationship 


J J Y na (0,A) Y kl (0,A) dcr 

cr 


r 4 tt if n = k and m = 1 
to if n^k or m/l 


( 2 ) 


with cr representing the unit sphere. The coefficients f nm are constants given by 

f«. = TZ rj F (©,X) Y nB (0,X) da j |m| ^ n, n^O (3) 

cr 

and constitute the spectrum of F . The error in the spectrum is thus 

e(fna) = JJ c (F(0,X)) Y n „(0,X) do ; \m\* H , n s0 (4) 

a 

There is at most a finite number of globally distributed measurements, which 
leads to the approximation 

e(f na ) M E c(F k ) Y na (0 k ,x k ) Act ; |m| ^ n, n > 0 (5) 

F k is the k-th measurement in the area element A a centered at ( 0 k , X k ) . It is 
assumed that all measurements are made in areas (blocks) of equal size and that 
the blocks uniformly cover the entire sphere (requiring polar satellites). Squaring 
the error yields 


e 3 (f n .) = tts E E MF k ) e(Fx) Y aa (0 k ,X k ) Y n# (0 1 ,X 1 )(Acr) 2 (6) 

16 TT k x 

The variance of the error is defined as the expectation of the squared error: 

m 2 (f n# ) = t^-sE E cov(F k ,Fi) Y a# (0 k ,A k ) Y^e^AO (ACT) 8 (7) 

16 TT k 1 

where the covariances cov ( F k , Fj) from the previous discussion are zero; i. e. , 
the noise is uncorrelated: 
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cov ( F k , Fi ) - 0 , k^l; cov( F k ,F k ) = m a ( F k ) 


(8) 


Then 

m 3 (f na ) 


— E m s (F k ) Y a 2 B (6 k ,X k ) (Aff) 3 
16 TT k 

m a ( F) E Act 


(9) 


if the variance of each measurement equals m s ( F) . Approximating the sum by 
an integral over the unit sphere and using (2), we finally arrive at the formula for 
the error spectrum 


m(f nB ) 


/ Ao“ 

“ V47T 


m(F); |m | £ n, n^O 


(10) 


Acr is written in units of square radians. 


We define the maximum resolution of the data by the frequency at which the 
signal to noise ratio is 1:1 . Hence, the spectrum of the signal must be known. 
Consider first the case of gradiometer measurements of the vertical gradient of 
the gravity disturbance, (9 2 T/^r 2 ) , where T is the disturbing potential and r 
is the geocentric radius ( r = R + h, h is the altitude of the satellite). In spherical 
harmonics, 

T(0,X,r) = Ry E (4T +1 t A nB Y„ B (0,X) (11) 

R is the mean radius of the earth, y = (kM/R 3 ) is an average value of gravity, 
kM is the product of the gravitational constant and the earth's mass, and the 
A nn are harmonic coefficients representing the spectrum of T/R y on the sphere 
of radius R . The zero- and first-degree harmonics are zero under the assumptions 
of identical centers, masses, and potentials of the level ellipsoid and geoid. Applying 
the radial derivatives, we find 


5 g T(8,X,r) 




_y 

R 


£5 /R \ n + 3 A _ _ 

E (A + l)(fr + J E A na Y n in( ® »^-) (12) 

n = 2 \ -k / e ~ — n 


Averaging the square of the n-th degree component of this signal over the sphere 
yields the degree variance (taking note of equation (2)): 




n ^ 2 


(13) 


The sum of squared coefficients A nm is modeled by Tscheming and Rapp (1974) as: 
£ A„ ffl 3 = - — - • - = 425,28 


b= -a 


y (n- 1) ’ 


2" i c n /„ _ 9 w„ (-999617) mgal , 


(n- 2)(n + 24) 
c 2 = 7. 5 mgal 2 


n > 2 


(14) 
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The c n are gravity anomaly degree variances. The degree and order variance 
of the gradient is then approximated by 


cr 


a 

na 



1 a /a*T\ 
2n + l n U 7 ') 


|ml s n, n ^ 2 


(15) 


The spectrum of the signal is consequently modeled as 

~ /3 S T\ _ 1 (n + l)(n+2) /R'\ n+3 /cT , . ^ 

MaTV " 72nTT J — R (t) ; lm| s “• " 4 2 (16) 

The degree n max for which the equality of o nB (5 2 T/Si^) and m( (P^T/Sr 3 ),!,,), 
given by (10), is attained then indicates the highest resolution in the data. 


The spectrum of the velocity can be determined by utilizing the law of con- 
servation of energy: 

I v 2 - V « E (17) 


V is the total gravitational potential (the minus sign is a matter of convention!. 


where 

v is the velocity, and E is the total (constant) energy. Decomposing V and 
into a systematic (reference) component and a signal part, we have 


i (v B + v) 8 - U - T = E 

i (v 8 2 + 2v,v + v 2 ) - U - T = E (18) 


where the reference velocity is defined by 


i v B 2 - U = E (19) 

U is, for example, the potential of a reference ellipsoid which at the satellite's 
altitude we approximate by the potential of a homogeneous sphere of mass M . 
Then from celestial mechanics we may show that 



Neglecting the second-order term, equation (18) becomes 

v = — T (21) 

v B 

The spectrum of the velocity signal, v , is then simply 

®'tiia( V ) — ® - nn( f ) j I Mil ^ n , n ^ 2 (22) 

where from (11) and (14) 

CT “< T > = 7In+T R fn^T) (7) S lml £ n, n * 2 (23) 
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We note that the corresponding degree variances are 

<7nV) = -4 <7n (T) ; n > 2 ( 24 ) 

CTb (T) = —^3 (■^•) n + 1 ; n * 2 (25) 

When comparing the velocity spectrum to the noise spectrum m( v nffl ) , 
equation (10), to determine the resolution of the data, it is assumed that the total 
velocity signal vector is measured. In practical situations, this ideal case is 
usually not realized, as for instance with satellite -to -satellite tracking in the 
high-low mode, only the component of the velocity signal vector along the line 
connecting the two satellites is detected. For a more thorough discussion of 
velocity measurements from satellite-to-satellite tracking, see e.g. Rummel 
(1978) or Hajela (1974). 

Now assume that two satellites are being tracked and the observations are 
the velocity differences between the satellites; for example, 4 v H = v P - v q , . 
where P,Q are the locations of the satellites travelling at the same altitude 
and on the same orbit in the low-low mode. However, we are now faced with 
a problem because, while the velocity v p , or v Q , under the approximations 
leading to equation (21) can be expanded in a spherical harmonic series, the 
difference Av H is a function depending also on a direction (azimuth of P 
with respect to Q) and on the central angle between P and Q and thus does 
not admit to an expansion of the form of equation (1). If we fix a direction and a 
separation, then it is possible to find, at least in theory, the spectrum of Av H 
(equ.(3)) despite its Singularities atthe poles (Hobson, 19 65, p.344). However, this 
spectrum is not merely proportional to the spectrum of the disturbing potential. A 
s imple r and not neces sarily less accurate course will be followed here. Cons ider the 
root meansquare (RMS) velocity difference over all possible direct ions 

A ^ H= fe/oJ (Vp “ v Q) 2da ] 2 (26) 

where a is the azimuth. This function is continuous over the sphere and for 
a fixed separation between P and Q, A v H can be expanded in a series of 
spherical harmonic functions. To obtain the power speetrumof Av H , we 
form the square of its RMS value over the sphere; 

cr 2 (AV h ) = 47 f J j Av h 2 dff 

. a 17 

= ~S J J (vp 2 + v q 2 - 2vp Vq ) dad CT 
= 2ct 3 ( v ) - 2cov(v P , vq) (27) 
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The last equation follows from the definition of the (physical) covariance func- 
tion (Moritz, 1972, p ,83) which, when the signal v is Isotropic and homo- 
geneous on the sphere, may be expanded as 


cov(v p , v'q ) = Z cr® (V) p a (cos^pq) 

a =3 

Substituting (24) into (28), equation (27) becomes (noting that cr s (v) 

00 

CT 8 (Av h ) - Z (l - P a (cos^rq)) cr*(T) 

a=S a 


(28) 

COV (V p , Vp)) 

(29) 


This is the sum of the degree variances of Av H , h e, its power spectrum 
components. The spectrum of Av H may be estimated from 


v a " ^Tl 4* 1 

„ /T 1 / 2(1 - P B (COS^po)) Cn 

Vy n - 1 V 2n+ 1 



m| ^ n, 
(30) 


n :> 2 


where (20) and (25) have been used. We then assume that the spectrum of 

the horizontal velocity difference in any one direction does not substantially 

differ from that of its RMS value over all directions; i. e. , cr a B ( Av H ) « or B B (Av H ) . 


The error of the velocity difference will be considered, not as the dif- 
ference of two errors, but as a directionally independent quantity to be rep- 
resented by the uncertainty of a single measurement, since the velocity dif- 
ference in satellite-to-satellite tracking is derived from a single Doppler 
measurement. Thus equation (10) is applicable if we assume that this uncer- 
tainty is a good approximation of the standard error of Av H , that is m( Tv H ) . 

None of the theoretical problems of the horizontal velocity difference 
occurs for the radial velocity difference; that is, when P and Q lie on the 
same radius. For the case of two orbiting satellites, such a radial velocity 
difference signal is a fictional quantity, however, since no two satellies at 
different altitudes could maintain an essentially constant separation. Never- 
theless, the analysis of the radial velocity difference may provide some con- 
firmation of that of the RMS horizontal difference if the signal is fairly 
homogeneous locally in space. 

The physical variance of the vertical velocity difference is 

ff 2 (Ar r ) = ^JJ(Av r fda 


= 2Cr 2 (v) - 2 cov (v P , v q ) 


(31) 
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by following the derivation of equation (27). The covariance between the 
disturbing potentials at two different altitudes is 

00 3 n + 1 

cov( Tp » Tq ) = E (tV) cr n s ( Tr ) P n (cos0pq) (32) 

n=3 \rp Tq/ 


where 0 „ s (Tr) is the degree variance of T on the sphere of radius R (taking 
r = R in equation (25)). Noting that at the altitude of point P , the degree 

( R 2 '? + ^ 

variance is cr n s ( T ) = (•— g-J cr n 3 (TR), and using equation (24), we find 


cov( vp , vq ) - ^ (-r) a n 8 (vp) Pa(cos(i)p,) , (33) 

Now 0 pq = 0 , since r P and r Q are colinear in the case of radial velocity differ- 
ence measurements. The subsequent substitution of (33) and (24), (25) into (31) 
then yields the degree variance of the radial velocity differences 


(Av r ) = 


(n-l)‘ 






(n- l) 5 


1 - “ 


&■) 1 

8 /rV +1 

{ o”) 

\rQ/ J 

\rp / 


(34) 


where r P < r Q . The spectrum of the signal may be estimated from 


O r nm(^V r ) /2n + 1 ^ ^ (35) 

Since the radial velocity difference, being independent of direction, is readily 
expanded in an harmonic series, equation (10) provides the corresponding error 
spectrum. 


The Quantities to Be Estimated 


The maximum resolution of the data, indicated by a signal-to-noise ratio of 
1:1 , is assumed to be also the maximum resolution of the subsequent estimated 
mean gravity anomaly or mean geoid undulation. Thus, the recovered mean 
anomaly or undulation contains no gravimetric information beyond the degree 
n max that corresponds to the maximum resolution. The error propagated into 
these estimates, due to the measurement noise, (i.e. the commission errort can 
be determined on the basis of the error spectrum model (equation (10) ). If we 
desire to estimate mean quantities, where "mean" denotes a simple average 
over a block (or spherical cap), then a further error, the truncation error, 
is introduced, since such a mean quantity by its definition comprises the entire 
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spectrum, although the averaging process dampens the magnitudes of the high 
frequency components. If, on the other hand, we desire to estimate a "filtered" 
quantity which by our definition contains no gravimetric information above degree 
n m ax» then obviously (by definition) no truncation error enters. The former case 
is adopted here. 


With the spherical approximation to the fundamental equation of geodesy 
(Heiskanen and Moritz, 1967, pp. 86-87) 

Ag = -T (36) 

Sr r ' 

the gravity anomaly, Ag , can be expressed in terms of spherical harmonics using 
(11) as follows: 

Ag(0,A,r) = y f (n-1) £ A na Y na (0,A) (37) 

n = 2 ' * ' m =s— n 

Also, from Bruns’ formula, N = T/y, the geo id undulation becomes (setting 
r = R in (11)): 

oo n 

N(e,X) = r E E A n ,Y nn( e,A) (38) 

n = 3 m =— n 

We define the mean gravity anomaly as the s imple average of Ag over a block on 
the sphere (which corresponds to current practice) and approximate this by a simple 
average over a spherical cap having the same area as the block: 


4g(9,X,B) = 2 T , ( fJoo S fe)JJ 4g dCT < 39 > 

where 0 O is the radius (spherical distance) of the cap O c and 2 tt( 1 - cos 0 O ) 
is its area. Substituting (37) (with r = R) into (39) yields the corresponding 
spherical harmonic expansion (see Meissl, 1971) 

oo n 

Ag(0,A) = y £ (n-l)jS n E A na Y nm (9,A) (40) 

n a: 2 o =— n 

Similarly, the mean geo id undulation is expanded as 

oo n 

N(0,A) = R E & E A na Y na (0,A) (41) 

n s 2 m n 

The coefficients j8 n effect a smoothing of the spectrum and are given recursively 
(Sjoberg, 1980) by 




S 3 


2 n-1 
n + 1 


cos 0o A,-i 


n - 2 
n + 1 




00 = It 01 = (1+cos 0o)/2 


(42) 
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From (12), the harmonic coefficient of degree n and order m of the 
vertical gradient of gravity at the satellite’s altitude is given by 



y_ 

R 


/ D \ n + 3 

(n + 1) (n + 2) A nB ; 


(ml £ n, n > 2 


(43) 


The corresponding coefficient of the mean gravity anomaly, from (40), is 

(^S)nm — y Ai(E — 1) -^nm 5 (44) 

and for the mean geo id undulation it is (see ( 41 ) ) 

(^)nm — R Aj -A nB \ (45) 

Substituting (43) into (44) and (45) , we obtain 


( )nm 

R ft, (n-1) / r \ n+3 /B 3 T\ 

(n + 1) (n + 2) \R/ Ur 3 / nB 5 

(46) 

(N)„ b = 

. R* ft /r\ n+3 /3 3 T\ . 

(47) 

y (n + 1) (n + 2) Vr/ Va^/ nB ’ 


The truncated mean anomaly and mean undulation are then 


Ag' 


aax 


R E 


n=3 



fl>(n-l) Y 

(n + 1) (n + 2) VR/ VBr s / nm nB 


(48) 



„3 n max n a ,_.n + 3 

IL y y & (jl\ ±^1 y 

y n =2 mi^n (n + 1) (n + 2) VR ) "" 


(49) 


TV 

Propagating the error in \y~p) na , given by (10), into these equations and averaging 
over the sphere (using (2)) yields the total commission error: 



(50) 

(51) 


Similarly, from (21) and (11), we find the (n , m) - harmonic coeffic ient of the 
velocity: 

( ^ )nm ~ R y -A nm i |m| ^ n, n s 2 (52) 

which when substituted into (44) and (45) gives 
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( v )na 5 


(53) 


( 4£v = v A. a . - -ii(X.y* 1 

(^g)nn v 8 R \R/ 


R \R 
— v„ $„ / r V* 1 

(N)„. - -“(*£•) (v) n . ; 

The corresponding commission errors, using (10) and (20), are 

n iuax ^ j 

“c( a 8) = [ n L^ 8(n " 1)2(2 n + 1) ('R 5r ) 

__ p /Act r naax /r 2 \ n+1 -i 

mc(N) = ^V47 m(V) [„| 3 S " a < 2n+1) ('F) ] 


,3 v n + 1 • 


(54) 


(55) 

(56) 


Recalling that <7 nB (Av H )» equation (30), is an estimate of the (n, m)- 
harmonic coefficient of Av H and that from (14), an estimate of 5 nm is 

A - - • H a n > ” * 2 < 57 > 

we obtain the approximation 

(AV H ) na = jyx JT( 1 - P n (cos^)p Q )) (— ) A M , |m| S n, n > 2 

(58) 


Substituting (58) into (44) and (45), it is not difficult to derive the commission 
errors as before: 


m 0 (Ag) = 
m 0 (N) = 



fyv V 4ir 


m(Av H ) 


m(Av H ) 



.^(n-..l) a (2n+.l) 

2(1 - P n (cos i/)pQ ) ) 

P 1 (2n + 1) 

2(1 -P n (CO8 0p P )) 



It is assumed in these equations that m(Av H ) m(Av H ). By not expanding 
Av h explicitly as a series of spherical harmonic functions, we were led to 
conjecture (equation (58)) that the spectral components of "Sv H are directly 
proportional to those of the disturbing potential. This unsubstantiated detail 
is not essential for the derivation of the commission errors, since we again 
average over the orders m. We require only that the corresponding degree 
variances are proportional to each other, this being the case as seen from 
equation (29). 


Finally, it is evident from (11) and(21) that the harmonic coefficient of 
the radial velocity difference is given by 


- 10 - 



<: n, n s 2 


(61) 


(Av r ) nm 




lm| 


where r Q > r P . Upon substituting (61) into (44) and (45), the total commission 
errors are readily found to be 


m c (Ag) 

mc(N) 





R 


rr— n aax 

^ m<4n) [.5,7 


a= 3 (1 -(V/r q f +1 ) 3 VtfV J 


(62) 

(63) 


The truncation errors, in each case, are simply (see equations (44) and (45)) 


fT (Ag) « 7 Y, & ( n_1 ) E 


n = n sax + 1 


(64) 


e T (N) = R E An E A„ b Y nB (65) 

n=n m «x+l ■ — n 

which when squared and averaged over the sphere become the RMS truncation 


errors: 


* 00 

m T (Ag) - [ E A^n] 

Ln=r >isax + 1 J 


m T ( N ) 


_ R 


— oo 

[2 

‘-n = n_ 


V +1 


- ftZ gy i * 

(n -1) J 


( 66 ) 

(67) 


where (2) and (14) have been utilized. 


A summary of the myriad of equations derived in this section is in order. 
The purpose was to obtain expressions for the errors in mean gravity anomalies 
and mean geoid undulations, the mean being defined as a simple average over a 
surface block. There are two types of errors; the commission error due to 
the uncertainties (noise) of the data, and the truncation error due to the finite 
resolution of the data (the mean by our definition contains all frequencies, 
though the higher ones are subdued). Equations (50) and (51) give the com- 
mission error if the data consist of gradio meter measurements; equations 
(55) and (56) correspond to the case of velocity measurements; if horizontal 
velocity differences constitute the data, then equations (59) and (60) describe 
the commission errors; while for the radial velocity difference, equations 
(62) and (63) are to be used. The truncation errors for each type of data are 
given by expressions (66) and (67). The total RMS errors are simply 

m T 0T ( 2 g) = [m 2 0 (Eg) + m ? r (Ag)p“ ( 68 ) 

m T oT(N) = [m s e (Nj + m 2 T (^]^ 
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The Data Density 


The parameter Act which enters the commission error is the least area 
element containing a single measurement. Considered as a block on the unit 
sphere, its dimensions are determined by the time interval between measurements 
along a track and by the (presumably uniform) track spacing on the equator after 
the duration of the mission. In the derivation of equation (10) it was assumed that 
the blocks Act are all equal in size and uniformly cover the entire sphere. Ob- 
viously, this assumption deviates from reality since the measurements increase 
in density with increasing latitude (for polar satellites). If h is the altitude of 
the satellite, then its period is approximately (Kepler’s third law) 



Given the duration of the mission, D , the number of orbits is 

D 

q = P (71) 

Assuming uniformly spaced orbital tracks after the duration of the mission, the 
across -track interval at the equator is 

TT 

di = q (radians) (72) 

noting that each orbit crosses the equator twice. The along-track interval is 
simply 

2 it P 

d s = s (radians) (73) 

where s is the sampling interval measured in units of time. The area element 
becomes (on the unit sphere) 

Act = d x • d a (radians) 3 ' (74) 

To obtain the intervals d x and d 2 in units of length, one simply scales them to 
refer to a particular sphere, e.g. on the earth, the intervals are Rdi and Rda. 
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Results 


A great number of tests can be conducted with the procedures outlined 
in the previous section. For initial considerations we examine the signal 
spectrum of the veloc ity difference and the vertical gravity gradient. The 
spectrum (signal /coefficient) is shown in Figure 1 (velocity) and Figure 2 
(gradient) for two different heights. In Figure 1 the horizontal and vertical 
velocity spectrum is shown for the case of the two satellites separated by 
300 km. Note that these two spectra differ by a factor of about 1. 4. 

Next we consider some specific GRAVSAT type missions. First consi- 
der a six month SST mission with a data spacing of 8 seconds and a data 
noise of ±lum/sec. For a radial separation of the two satellites of 300 km 
the anomaly and undulation accuracies to be expected are shown in Table 1 
for two different altitudes. In this table accuracy estimates from both a hor- 
izontal and a radial separation mission are given. It is clear that there is 
not a significant difference in the results. 

Similar results are seen from a g radiometer miss ion as shown in Table 2. 

The results presented in both Table 1 and 2 are based on the Tscheming/ 
Rapp anomaly degree variance model. Computation made with the more 
recent two component model of Rapp (1979) yield values about 20% larger than 
what is given in the tables. 


Table 1 

Accuracy of Mean Anomaly and Mean 
Geo id Undulations as Determined from 
SST Tracking (Noise = iljum/sec, | Radial Mission ; 

( Horizontal Mission 


Separation = 300 km) ; T/R, 208 Degree Variance Model 



H = 150 km 

H = 180 

km 

Block 

Size 

Anomaly 

■ 

1 


Undulation 

2° 

0 . 9 mgal 

1.6cm 

1,2 mgal 

2. 6cm 


0.8 

1.4 

1.1 

2.4 

1° 

2. 2 mgal 

4. 1 cm 

3. 3 mgal 

7.3 cm 


2.0 

3.6 

3.0 

6.3 

30' 

9. 1 mgal 

16. 2cm 

10. 6 mgal 

21. 1cm 


8.7 

15.2 

10.3 

19.8 
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Table 2 

Accuracy of Mean Anomaly and Mean 
Geoid Undulations as Determined 
from Radial Component Gradiometry 


(Noise = 10.01 EU.) 



H = 150 km 

i 

H = 180 km 

Block 

Size. 

Anomaly 

Undulation 

Anomaly 

f 

Undulation 

2° 

0. 8 mgal 

1.9cm 

1.2 mgal 

2.9 cm 

1° 

1.9 mgal 

3.7 cm 

3.1 mgal 

7.2 cm 

30’ 

8. 5 mgal 

14.5 cm 

10 . 4 mgal 

20. 6cm 


Table 3 

Truncation and Commission Error 
for H = 180 km and 1° Blocks 




Commission 

Error 

Truncation 

Error 

Total 

Error 

Type of 
Measurement 

Njjjax 

Ag N 

Ig N 

Ag W 

Gradiometer 

275 

1.9 mgal 5.4 cm 

2. 5 mgal 4.7cm 

3.1 mgal 7.2cm 

Radial 

Velocity Diff. 

266 

1.8 mgal 4.8 cm 

2. 7 mgal 5.5cm 

3. 3 mgal 7.3 cm 

Horizontal 
Velocity Diff. 

273 

1.3 mgal 3.9 cm 

2. 5 mgal 4.9 cm 

3.0 mgal 6.3 cm 


Additional tests have also been made with a mission having the satellites 
at an altitude of 160 km with various range rate accuracies specified over 
a 4 second integration time. Some representative numerical results are 
given in Table 4 for the accuracy of die 1° x 1° anomaly and undulation 
recovery. -16- 





















Table 4 

1° x 1° Anomaly and Undulation Accuracies 
for a Low-Low 6 Month Mission 
(Horizontal Separation = 300 km), Height = 160km 


Noise 

m(Anomaly) 

m( Undulation) 

1 jim/sec 

t 2. 1 mgal 

- 3. 7 cm 

2 

2.6 

5.2 

4 

3.2 

7.0 

5 

3.4 

7.7 

6 

3.6 

8.4 

8 

3.9 

9.5 

10 

4.2 

10.4 

50 

6.9 

21.1 


Comparing Table 1 and Table 2 we see very similar results. We also 
give in Table 3 a comparison of the commission error and the truncation 
error for the radial and horizontal measurement and the radial gradiometer 
measurement for H = 180 km. 


From the table we see that the estimated accuracies are not linearly depen- 
dent on the data noise. This was also seen in the results of Rapp and Hajela 
(1979) using the method of least squares collocation for anomaly and undulation 
accuracy estimates. 

Computations were also made to obtain results to compare to that given 
for a high/low case in Table 3 of Rapp and Hajela (1979). The accuracies 
obtained from the method described in this paper were approximately 20 % 
smaller than those found through the least squares collocation techniques 
used by Rapp and Hajela (1979). This is not unreasonable as the collocation 
results used less data than implied by the techniques of this paper. 

In these results we have concentrated on looking at accuracy estimates 
in blocks. It is also possible to look at the accuracy estimates by degree n. 

To do this we can look at the percentage error up to the cut off degree. This 
is done by dividing the accuracy at a given degree by the expected signal. This 
relative accuracy estimate will hold for any quantity of interest (potential 
coefficients, anomalies, undulations etc.). 

Such a computation has been made for a low -low 6 month radial mission 
with H = 160km, and a 4 sec integration time. Results are presented in 
figure 5 for two noise estimates l£lm/sec and lOjUm/sec. These percentages 
can also be converted to undulation and anomaly errors. Using the Tscheming/ 
Rapp anomaly degree variance model, and the Ijum/sec noise level of the 
mission in the above paragraph, the undulation error at degree 300 was 2. 1cm 
and the corresponding anomaly error was 1. 0 mgal. 
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The percentage errors by degree were also computed for the low-low 
horizontal mission. These results are also plotted in Figure 3 where we 
see only a small shift to the higher spherical harmonic degree. 

One additional way to look at the results is to consider the cumulative 
commission error out to the maximum degree that can be determined. This 
procedure would correspond to the accuracy of the computation of a point 
quantity with a truncated series. Some results for the 160 km, 6 month low- 
low horizontal mission are shown in Table 5 using the Tscheming/Rapp 
anomaly degree variance model and a noise of 1 Mm/sec. 


Table 5 

RMS Error of "Point" Computation 
Assuming No Truncation Error 


N Ba x 

m{ Anomaly) 

m( Undulation) 

50 

t . OQOmgals 

t ,007cm 

100 

.003 

.022 

150 

.026 

.123 

200 

.128 

.460 

250 

.584 

1.636 

300 

3.195 

7.327 


Summary and Conclusions 

This report has examined the accuracy in the determination of the earth’s 
gravity field, to be expected from low-low range rate or a gradiometer satel- 
lite mission. The theory used is one of error propagation and truncation effect. 
We started w ith the type of measurement and computed the noise effect of that 
measurement on the spectrum of the earth’s gravity field. The maximum de- 
gree that could be determined was set when the propagated noise equaled the 
estimated signal based on some model of the spectrum. We then examined the 
computation of mean anomalies and mean geo id undulations by considering two 
error components. The first is the propagated commission error due to the 
noise in the measurements and the second is the effect of the neglected higher 
degree terms on the mean anomaly and mean undulation. In doing this we used 
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the Pellinen smoothing operator to relate a spherical harmonic expansion to 
a mean value. Such a procedure closely approximates that actual procedure 
used to determine a mean quantity in a Hock. 

It should be noted that the procedure developed in this report is not an 
optimum estimation procedure such as is used by Breakwell. However if 
our results are compared with those obtained by Breakwell using the same 
filter or weighting parameter, they are within 15% showing that our procedure 
and the optimum estimation procedure do not yield significantly different results. 

Our procedure basically ignores earth rotation and simply assumes that 
data will be obtained on a global basis on some interval dependent on the lifetime 
of the mission. In addition we have assumed no orbit error is present in the 
data. 


In carrying out these developments we assume that the spherical har- 
monic expansions are valid everywhere on and above the geoid. In fact this 
is not the case because of the topography and so some additional attention 
needs to be placed on the downward continuation problem. However we have 
not discussed in this report a data reduction procedure but an error determin- 
ation procedure. When the actual data reduction is carried out the downward 
continuation question must be considered if results are to be given on the 
earth's surface or the geoid. 

The equations derived were applied to several test missions. We found 
that a low-low mission with the two satellites having a horizontal separation of 
300 km at a height of 150 km and a data noise of * lum/sec gave essentially the 
same results as a radial component gradiometer (±0. 01E) at the same altitude. 

Tests were also done with a radial low -low 6 month mission, H = 160km, 
with various data noise. For ljUm/sec noise (4 second integration time), 
the 1° x 1° anomaly would be determined to ±2. 3 mgal while the corresponding 
undulation would be ± 4. 3 cm. If the noise was raised to 6 /hn/sec the cor- 
responding values would be ± 4. 1 mgals and 10.0 cm. 

The results presented in this report clearly show that significant im- 
provement in the determination of the earth's gravity field can be expected 
if a gravity mapping satellite mission could be flown. The equations developed 
here have been programmed in a very fast error analysis study program that 
can be used for many other test cases. We now need to develop a suitable 
data reduction procedure. 
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